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ABSTRACT 


The Monterey-Miami Parabolic Equation (MMPE) Model is a full-wave 
underwater acoustic propagation model that utilizes the split-step Founer marching 
algorithm. Previously the MMPE model was implemented in Fortran language and ran 
with a simple command line interface either in a Unix or DOS command window. After 
the Fortran code was run, the resulting binary data output file was post-processed using 
Matlab routines to extract specific field data and present the results in graphical form. 
This approach requires the user to have installed both Matlab and Fortran compilers. The 
MMPE model and associated acoustic processing tools are now rewritten in the object- 
oriented language Java. This new version of the MMPE model built within a Windows 
framework is called WinMMPE. Integrating the model, the post-processing calculations 
and the graphics generation together with a graphic user interface has produced a more 
attractive tool for users. A user-friendly, efficient, and accurate full-wave acoustic 
propagation model with enhanced visualization can make it easier to assess the spatial 


transmission loss in the underwater acoustic environment. 
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I. INTRODUCTION 


The parabolic equation method for predicting underwater acoustic propagation 
was introduced to the underwater acoustics community by Tappert (1977). The parabolic 
equation (PE) is a one-way wave equation approximation to the full Helmholtz wave 
equation that defines the continuous-wave (CW), single frequency acoustic field in 
general media. Because the PE is still a full-wave equation, it includes all of the full- 
wave effects (e.g., diffraction) that are neglected in more rudimentary ray-based models. 
It is also well suited to handle range-dependent media as it naturally computes the 
coupling between vertical normal modes of the waveguide. Additionally, the one-way 
formulation of the PE allows for implementation with highly efficient numerical 
algorithms. The most common methods are the split-step Fourier (PE/SSF) method 
(Hardin and Tappert, 1973), the implicit finite difference (IFD-PE) method (Lee, et al., 
1981), and the finite element (FEPE) method (Collins, 1990). 


Of these methods, the split-step Fourier method was adopted in developing the 
University of Miami Parabolic Equation (UMPE) Model. The UMPE Model was 
implemented by Smith and Tappert (1993). Since Prof. Smith moved to Monterey in 
1995, the UMPE Model was upgraded to incorporate some improved numerical 
techniques and rewntten. This new version was then named the Monterey-Miami 


Parabolic Equation (MMPE) Model. 


The MMPE Model was written in Fortran and implemented using a command line 
structure available in either a Unix or DOS command window. A Fortran compiler is 
required to create the executable program. Its output file is a binary data file composed 
of the complex values of the acoustic field function for various grid points in range, 
depth, and azimuth for multiple frequencies, in general. From this data, acoustic pressure 
and transmission loss can be extracted for a variety of simulated receiver locations within 
the computational grid. This post-processing of the binary data file and the subsequent 


graphical manipulation of the results are performed using Matlab. 
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Many underwater acoustic models are wnitten in Fortran language. The compiled 
executable files run and output data files of the acoustic prediction are produced. This is 
then processed using a graphics package to display the results of interest. Such command 
line interface programs have little advantage in Windows operating systems. First of all, 
users have to handle the programs separately. Most users would probably prefer using an 
integrated software package because it is easier to handle. Second, the program’s 
command line interface and the Matlab command window are text based. Users are 
required to type text in the command line. It is not a user-friendly computing 
environment. Third, in standard Fortran the size of data arrays must be declared before 
they are used. Developers must always consider the maximum array size that may be 
needed. This can sometimes make the programs waste unnecessary memory depending 
on some input data, or users might want to use more input data than the maximum 
declared. Fourth, the Fortran program is a stand-alone program. It is implemented in an 
independent computing environment. It is not always easy to combine with other 
software or systems, and it requires a Fortran compiler (not standard on all computer 


systems). 


To address these issues, the MMPE model must be re-programmed. It should be 
an object-oniented program, and the propagation model and post-processing graphical 
program should be integrated as one. It should have a new graphic user interface (GUI) 
and should be convertible to any kind of system. For these reasons, the choice for the 
programming language is Java. Java is an object-oriented language and network 
language. Java can dynamically create the data array size because it treats a run time 
object system. It is also possible to be multi-threaded, which will allow it to run several 
tasks in parallel within the same program. Lately, the Java language has become very 
popular among many programmers since it was announced in 1995. It 1s open-source and 
freely served by Sun Micro Systems Corporation. More and more class libraries have 


been and continue to be developed. 


The main goal of this thesis is to develop a user-fnendly, efficient, and accurate 


full-wave acoustic propagation model with enhanced visualization. This will be 
Z 


accomplished by integrating the MMPE Model and graphic functions into a single 
program with a new graphic user interface (GUI). The integrated software will make it 
easier for users, such as the war fighter in an operational situation, to assess the acoustic 


environment of the battle space. 


Chapter II provides a summary of the theory of the MMPE model and its 
implementation and graphic manipulation. It derives the general form of the parabolic 
equation to be used and describes the algorithm used in the MMPE model. The post- 
processing and graphic manipulation routines are also defined. In Chapter III, the design 
of the Java-based WinMMPE and its functions are presented. The details of the GUI are 
also presented. The model implementation, post-processing, and graphic manipulation are 
described. In Chapter IV, the output results of the WinMMPE are described. Sample 
outputs of WinMMPE are compared with data from the MMPE Model. This chapter also 
presents a discussion of the efficiency in terms of main memory used for implementation 
and running time. Finally, Chapter V provides general conclusions and recommendations 
for further work. Suggestions for methods of improving efficiency and the addition of 


new modules are provided. 
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Il. MONTEREY-MIAMI PARABOLIC EQUATION (MMPE) MODEL 
A. INTRODUCTION 


This chapter provides a general description of the Monterey-Miami Parabolic 
Equation (MMPE) Model, its algorithm, and its implementation. All contents in this 
chapter are summarized from the report of the University of Miami Parabolic Equation 
(UMPE) Model by Smith and Tappert (1993) and the recent article of Smith (2000) that 
analyzed the accuracy of the MMPE Model. Accordingly, the nomenclature used in this 


chapter is the same as in those reports. For more details, the reader should refer to those 


papers. 


B. THE GENERAL THEORY 


We begin by considering a general time-harmonic, continuous-wave (CW) 


acoustic field at angular frequency Win a cylindrical coordinate system (r, Z, @), 

PAC aZn ORG) = 9 (nize eam (2.1) 
Combining Eq. (2.1) with the inhomogeneous wave equation for a point source (Jensen, 
et al, 2000), 


2 
Dy, ‘aie ee = -476 (Xx -X,) , (2.2) 
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leads to the Helmholtz equation for the CW acoustic field, 





iG) sen i ledaneecua ae 
a Rear ee eeou Kot (.2-Cu4P Olea— x.) . (2.3) 


In this equation, k, is the reference wavenumber defined as 


4y) 
k,=—., 
0 a (2.4) 


C, 18 the reference sound speed, n(r,z,@) is the acoustic index of reference defined as 


Co 
n(r, Z,Q) = ezP) ; (2:5) 


and c(r,z,@) is the acoustic sound speed that vanes throughout space, in general. The 


source function is defined as a point source with reference source level P, at a reference 


distance of R, =1m, and 


§(z) = —d(e-z, \5(r) , (2.6) 


where r = 0,2 = 2, refer to the position of the source. Solutions to Eq. (2.1) may 


be computed over a bandwidth of frequencies, and the results can then be Fourier 


synthesized to produce predictions of the arrival time structure due to broadband pulse 


propagation. 


In most ocean environments, the influence of the azimuthal coupling term is 
negligible and may be ignored. To further simplify Eq. (2.3), we shall also neglect the 


source term initially. We may then define 


[ey = sul ; (2.7) 


and substitution into Eq. (2.3) yields the free-field wave equation 
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=e. 
Or* dz’ -} (2.8) 


Next, we introduce the operators 





) 
Po=— 
~~ (2.9) 
and 
1/2 
lead 
=|n°>+———| . 20 
Ox, . ke a ( ) 
In the far field (kyr >> 1), Eq. (2.8) then becomes 
(P2 +420? W=0. (2.11) 
Expanding Eq. (2.11) into incoming and outgoing wave equations, we obtain 
Ou 
ees | = 
ee 5, Lon (2.12) 


for the outgoing wave when the incoming wave is negligible. Now u can be defined in 


terms of the PE field function, w(r,z), as 
u=yw(r,ze . (2.13) 
Substituting this into Eq. (2.12), we finally obtain the PE for the field function, 
a aa = -ik,H 
= oV nick 2 pV ae 0 a ’ (2.14) 
where 


ao (2.15) 


is an Operator having Hamiltonian form. This operator defines the evolution of the PE 


field function in range. 


C. THE SPLIT-STEP FOURIER (SSF) ALGORITHM 


We introduce the PE field function, y, by formally defining the relationship with 


the acoustic pressure field as 
Ro -1/2 ikor 


normalized such that at r= R), lvl =] and | p| = P,. The computational evolution of the 
field is now accomplished via a marching algorithm of the form 


y (r+ Ar)=P(r)y(r) . (2.17) 


From Eq. (2.14), the propagator ®(r) can be shown to have the approximate form 


@(r) = go tkoHop (rar (2.18) 
where 
—_ l ert ; ? 
H,,()=— © dr'H,,(r’) . (2.19) 
if r 


If the range step Ar is small enough, the operator may be assumed to be relatively 


constant over a range step. Evaluating the operator at mid-step, for instance, gives 


H,,(r) = Hal +5 (2.20) 


From Eq. (2.10), it is observed that the operator Q,,, and therefore Ho,, are 
pseudo-differential operators due to the square root. An approximation must then be 
made for numerical implementation. In order to implement the resulting parabolic 
equation using the split-step Fourier method (Hardin and Tappert, 1973), it is necessary 


to separate the index of refraction term from the vertical derivative term, 1.e., 


eet UT |. (2.21) 


It has been shown (Jensen, et al., 2000) that a solution which is second order accurate in 


range step size, Ar, is obtained when the propagator function is defined as 


Rees npg 
=k ——U op (r+Ar) —iky ArT, Ike op (r) 
é 


O(r)=e 2 e , (2.22) 


Using the wide-angle operator splitting introduced by Thomson and Chapman (1983), we 


define 





Tr. 21-\1¢2 2.23 
op le dz? ( : ) 


and 


U,, =-(n-1). (2.24) 


Because U,, is a scalar operator, it forms a diagonal matrix which can simply be 
multiplied by the PE field function at each range step. However, 7,, is a differential 


operator, is not diagonal, and couples the solution between depth points. The split-step 


Founer algorithm overcomes this by recognizing that the corresponding operator in the 


vertical wavenumber domain, 7,,, 1s diagonal and can be treated using a simple 


multiplication. 


In the PE/SSF implementation, generally, a multiplication of the z-space operator 


Ar 
Tio YaplT) apnea 
"2°""" is defined at the initial range-step. The result is then transformed to the 


vertical wavenumber, k,, domain. After a multiplication by the k,-space operator 


-ikgArT : ‘ : 5a . 
e "the solution is transformed back to the z-domain where a multiplication of the z- 


Ar 
~iky—U (r+ Ar) ; 
space operator e 7? ” is defined at the end of the range step. A fast Fourier 


transform (FFT) algorithm is used to accomplish the transformation which assumes a 


convention 


y(z) = FFT(W(k,)) . (2.25) 
The representation of the PE/SSF implementation is then 


A 
~iky——U op (r+r.2) Se UE 


w(rtAr,zj=e 7° cFer| en ee vino) (2.26) 


where, in k_-space, 


5 WV? 
: k. 
ji =-}-{) | | (2.27) 


D. THE BOUNDARY CONDITIONS 


The ocean surface in the MMPE model is defined by a Dirichlet boundary 


condition consistent with a pressure release boundary, 


Vie '0) =O (2.28) 


For this boundary condition, an image ocean method is used whereby an exact copy of 
the environment is defined for negative depths and the field function for these depths 1s 


defined by 


W(-z) = —w(z) . (2.29) 


Accordingly, the data array 1s twice as large as that needed to describe the real acoustic 


field. 


In the MMPE model, shear wave propagation in the bottom is treated by an 
equivalent fluid technique. Thus, the propagation of shear waves are not computed but 
the effect of shear on the boundary reflection coefficient is included. Specifically, the 
equivalent fluid approximation of Tindle and Zhang (1992) is employed at each bottom 


interface. 
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It is assumed that the interface between the bottom of the water column and the 
top of the basement, or sediment, layer is characterized by a discontinuity due to a sharp 


contrast in sound speed. Specifically, the sound speed at the interface is defined by 


e(z) =c,,(z) +[e,(z)—¢, (DIN (z-z,) ; (2.30) 


where c, and c, are the sound speed of water and bottom, respectively, and the 


Heaviside step function is defined by 


0, ¢€<0 
H(g)=4Y2, ¢=0 , (2.31) 
| ore Ge (6) 


where ¢ = z-—z,. In order to implement within the numerical scheme, it is necessary to 


approximate this by a smooth, continuous function. For that purpose, this function is 


approximated by a hyperbolic tangent function, 


H(f)= 3p 4 ant = | (2.32) 


c 





or, equivalently, 


H(f)=(l+e%/*) , (2.33) 
where L, is the sound speed mixing length. 


In order to create the most realistic interface condition for reflections from a 
sound speed discontinuity, Smith (2000) showed empirically that the most accurate 
results were obtained by defining L, as a fraction of a wavelength, specifically A/10. For 
sampling purposes, the default minimum value for L, is set at L..;,,=Az. Thus, 


L. = max(Az,A/10), so the most accurate results are obtained efficiently when 


Az < A/10. 


1B 


We must now consider the effect of the bottom density contrast, something which 
was neglected in the definition of the wave equation given in Eq. (2.2). A proper 
treatment including density produces the same expressions given above if we replace the 


index of refraction n by an “effective” index of refraction, 


SO | ee Cla 
n° =n +—~|—V°p--|-—V , 
ve E p | e) (2.34) 





Then, the pressure field p(r,z) is defined by 


R - ikgr 
p(r.2)= P| “Ommuaae (2.35) 
Por 


where , is a reference density defined by p, =1.0 g/em®. The environmental 


propagator function in Eq. (2.24) must now be defined as 


OF) CAG rez) (2.36) 


where U,(z) is the same as the previously defined environmental potential function in 
Eq. (2.24) and U.,(z) is a new operator that includes the effect of the density 


discontinuity. 
If we assume that 9 = 9(z) only, the density is defined by 


DQ) 10, Fs = (0), Nl eee (2.37) 


where 9,, and 9, are the density of water and bottom, respectively. As before, we must 


replace the Heaviside step function by a smooth, continuous function scaled by a density 


mixing length, L,. Within the MMPE Model, the function U2 1s approximated by 
Sg p PP 


UI oo aye (2.38) 


U3) =- : 
° is aeree 


where 
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_ 1-(p,/p,)” 
1+(0,,/0,) 
For the density mixing function, a cubic spline over the finite interval 


-L,<¢SL, is defined. Then H(¢) is defined by 


0 gor 
ji, 
= noe | 5 Se 
3 L, s 2 
we eG a Gp uf 
H(¢)=,—+—-- , oe << -— CC. 2.40 
Lf #4 sce. 
ie 
peep o , —<OSL, 
3 jb 2 
i sy 
The first derivative of this function is 
ta by. 
oo 4) , ee C<-—+ 
=, L, Z 
a . 1 C L, L, 
A =0(¢) =; —11-4 => ge Meee ee 2.41 
)= 3) =} 7 {£] tall (2.41) 
DE ra -£) a 22 ib, 
L a: 
0 ; Bel, 
where 0(¢) is the smooth approximation to the Dirac delta function. The second 


derivative is then 
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Equations (2.38) and (2.42) provide the necessary expressions for computing the density 


potential function with a cubic spline polynomial smoothing function. 


In order to create the most realistic interface condition for reflections from a 
density discontinuity, Smith (2000) showed empirically that the most accurate results 


were obtained by defining L, ~ 2A. For sampling purposes, the default minimum value 
for L, is set at L, 4, =5AzZ. Thus, L, = max(5Az,2A), so the most accurate results are 


obtained efficiently when Az < 2A/5. 


E. THE SOURCE FUNCTION 


Near the source, we may assume that the effect of density is negligible and the 


operator in Eq. (2.16) is near unity. The acoustic pressure near the source is then defined 


| R. . 
p= ee (2.43) 


where KR, is a reference range defined by 1m and P, is the source pressure at R,. The 


simply by 


source level, SL, becomes 


r 


P 
Si = 20108 | dBreP.,R,. (2.44) 
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where P. is a reference pressure, typically defined by P. =1 wPa. 


r 


From Eq. (2.43), w(r =0, z) can be expressed as 


|e 
= 0,z) = lim— —- RZ) i. 
w(r Zz) OP VR, P(r, Z) (2.45) 


In the vicinity of a point source, the pressure field takes the form of a spherical Green’s 


function, 1.e. 


p=——e® | R=avr?+z2? , a=PR,, (2.46) 


P. 
where a is defined by requiring | D| = = at R=R,. If we then define the PE field 
1 


Starting condition for a point source at depth z,; by 


w(r =0,z) =@0(z-Zs) ; (2.47) 
then it can be shown that the proper normalization is 


a= (2.48) 





To represent this field in the vertical wavenumber, k,-domain, we perform a Fourier 


transformation of Eq. (2.47) to produce 


Wir =0;k)=—Ziasin(k.z,) (2.49) 
which also includes the influence of the image source. 


For the wide angle approximation in the vertical wavenumber domain, we need to 
add two factors in Eq. (2.49). First of all, to produce the correct solution in the far-field, 


-/4 
k Z 

the factor [ = - (k.| <k,) needs to be added to the wide angle source (Thomson 
0 
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and Bohun, 1988). Secondly, in order to take advantage of the MMPE depth grid 


convention, defined by 


(2.50) 


pies 
we need to add a phase term in the wavenumber domain of e 7. Finally then, the PE 


starting field in the k.-domain for the wide angle source is given by 








0 0 


W(r=0,k.) =—2i “to sin(k,z.} 1 : = (2.51) 
r= : o =-—-Z i — - e - , i 
: 27k ~~ k- 


F. IMPLEMENTATION 


There are two aspects of the implementation of the MMPE Model. The first is the 
propagation model implementation that produces a binary output file. The second aspect 
is then the post-processing of this binary data file and the graphical representation of the 
results. This is treated using Matlab routines. The results of the output PE solution will 
be discussed in Chapter IV. In this section, the run-time process and the output structure 


is discussed. 
1. Model Implementation 


The MMPE Model can be run using a DOS command line interface as illustrated 
in Figure 2.1 (a). Prior to execution, several ASCII input files must be defined. The 
primary input data file is named “pefiles.inp’. Within this file, all other filenames that 
provide details of the environment (sound speed profiles, bathymetry, etc.) are defined. 
Of all file names, only the specific name “pefiles.inp” is required. Also within this main 


file are specifications for the output data file (number and range of output depths and 
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ranges) as well as optional declarations of computational parameters (FFT size, range and 


depth grid sizes, etc.). 


When run, the propagation model asks the user to choose either an accurate mode 
or an efficient mode of calculation (unless the grid sizes have been explicitly specified in 
“pefiles.inp’’). In the former, the depth grid is 1/10 of a wavelength and the range erid is 
a wavelength. In the latter, the depth grid is half a wavelength and the range grid is three 
times a wavelength. The model calculates the PE field depending on the user’s choice. 
During the calculation, the range step displays on the screen, along with the current 
frequency and radial bearing being computed, as illustrated in Fig. 2.1 (b). The result of 


the calculation is the binary output with filename defined in “pefiles.inp.” 


The structure of the output binary file is shown in Figure 2.2. For all but the final 
records (containing bathymetry information), each record length is twice the number of 
points stored in depth. The factor of two accounts for the real and imaginary parts of the 
complex PE field function. The first record contains header information providing details 
about the remainder of the data and the calculation as follows: the total length of the data 
file; the reference sound speed; the number of frequencies computed; the center 
frequency; the frequency bandwidth; the number of range points output; the minimum 
range for output; the maximum range for output; the range step size used in the 
calculation; the number of depth points output; the minimum depth for output; the 
maximum depth for output; the number of radials; the angular difference between radials; 
the bottom depth at the source location; the deep bottom depth at the source location; and 


the source depth. 


For multiple frequency and multiple radial calculations, the model first computes 
the solution out in range for the initial frequency and radial. This begins by computing 
the starting field, which is the same for all radials at that frequency. It then computes the 
solution for all subsequent radials at that frequency. After computing the solution out in 
range for every radial, it then begins the calculations for the subsequent frequencies. The 


result of the ordering is the binary data file organization depicted in Figure 2.2. 
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Figure 2.1: The MMPE model run-time environment: (a) initiation of model run; 
(b) model output while running: 
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Figure 2.2: The structure of the MMPE output binary data file. 


Note that the data ordering places adjacent radials next to each other rather than 
adjacent range steps as might be expected from the actual run-time ordering. This was 
done to maintain consistency with other versions of MMPE, which incorporate true 3D 
effects with azimuthal coupling. In those versions, all radials are computed at each range 
step. Also note that the complex PE field function values are stored as real and 


imaginary data pairs. 
2. Post-Processing and Graphic Manipulation 


The post-processing implementation runs in the Matlab command window. The 
initialization routine is called “peout1”, which requests the user to enter the name of the 
binary output file of the MMPE Model to be processed. After entering the file name, the 
data of the output file is loaded to the variables in “peout1” as shown in Figure 2.3 (a). 
The main processing routine, “peout2”, can then be run multiple times without re-running 


the initialization routine so long as the same file is being processed. After typing 
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“peout2”, six options for processing are given, as shown in Figure 2.3 (b). The user then 
chooses the option for the desired processing and data extraction. For the option chosen, 
all acoustic quantities of interest are computed, and a corresponding transmission loss 


field is plotted. The outputs will be discussed in more detail] in Chapter IV. 
For all output options of interest here, transmission loss (TL) is defined by 


|p(r,z,9)| 
P 


0 


TL(r, Z,Q) = 01s | dB re lm, (2.52) 
where Pp is the source pressure amplitude measured at the reference distance of R, =1m. 


From the definition of the pressure field in terms of the PE field function in Eq. (2.16), 


this may be rewritten 


] 
TL(r, z,Q) = 20log,,| —=lw(’, z,Q) 
43 | | (2.53) 


= 20log,.|w(r, z,9)| -10log,,r dBrelm 


Thus, the solution of the PE field function provides all the information necessary to 


compute TL with the exception of the simple cylindrical spreading factor. 
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Figure 2.3: Post-processing and graphic implementation: (a) loading the MMPE output 
file; (b) choosing one of the processing options. 
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Ill. WINDOWS MMPE (WINMMPE) 


A. INTRODUCTION 


As described in the previous chapter, the MMPE model can be executed via Unix 
shell or DOS command. Then, the post-processing and the graphic manipulation are 
implemented in MATLAB. Since the final goal of this thesis is integrating those 
functions and providing a better interface for users, a Windows version of the MMPE 
sausfying these goals has been created. It is called WinMMPE. This chapter describes 
its design and functions. Also, its implementation 1s described in terms of the model and 
the graphic manipulation. The WinMMPE application is executed by clicking its icon 


(Figuve sar). 





Figure 3.1 — Icon of WinMMPE 


B. DESIGN AND FUNCTIONS 


WinMMPE is composed of five frames, consisting of the main frame, model 
contro] frame, input data editor frame, graphic control frame, and graphic output frame. 
Each of the frames has its own design and functions. This section describes how they 


have been designed and what kinds of functions they have. 


1. Main Frame 


A double click of the left mouse button on the icon makes the WinMMPE display 
the main frame (Figure 3.2). It provides a full Windows graphic user interface for all 
WinMMPE functionalitv. There are five menus, which are “File” menu, “Edit” menu, 
“View” menu, “Window” menu, and “Help” menu. Each menu has several menu items 


and all are shown in Table 3.1. In the “File” menu, there are five menu items. “New 
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Source” menu item renders model contro] frame with main input file, “default.inp.” 
“Open Source” menu item renders an open file chooser. Users can choose a “.inp” file 
as a main input file. If other formats are chosen instead of “.inp,” file errors will occur. 
“Save” and “Save As” menus are used to save the main input file with the given file 
name, or to save it with another name if desired. “Exit” menu item makes WinMMPE 
terminate. The “Edit” menu has three menu items: “Cut,” “Copy,” and “Paste.”” When 
editing input files in the “Input File Editor” frame, these menu items are supposed to 


work for the frame. However, they do not work yet and are left for future work. 








FS windows Dorabelc need Cee 


File Edit View Window Help 









Figure 3.2 - WinMMPE Main Frame 
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New Source 








Trew sone 


Cut 


Pe a ee Frame 
View eee File Editor Frame 


| Graphic Control Frame —_| Control Frame 
Opened Frames 


Web Help 


About Wi 
Table 3.1 —- Menu and Menu [tems of Main Frame 





Help 









In the “View” menu, there are three menu items: “Model Control” frame, “Input 
File Editor” frame, and “Graphic Control” frame. Whenever any of these menu items is 
selected, a file chooser is opened. For “Model Control” frame and “Input File Editor” 
frame, it shows “.inp” file as option. For “Graphic Control” frame, it shows “.bin” file as 
option. The “Window” menu shows the opened frames as menu items. Users can check 
how many window frames are on the main frame and can pop up the frame they select. 
The “Help” menu has two menu items: “Web Help” and “About WinMMPE.” The 
former lets users give the information about menu and each menu item of WinMMPE. 


The latter is the general explanation of WinMMPE. 


The nine buttons in the tool bar below the menu bar are seen in Figure 3.1. From 
left to right, they are: “New Source” button, “Open Source” button, “Save” button, “Cut” 
button, “Copy” button, “Paste” button, “Model Control” frame button, “Input Data 
Editor” frame button, and “Graphic Control” frame button. Each button is able to work 


like the corresponding menu item previously described. 
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2. Model Control Frame 


This frame launches the main input data file and the sound source input file, 
which have the same form as the ones used in the MMPE model. It comprises seven 
panels: “Output Data Files” panel, “Depth Info” panel, “Range Info” panel, “Basic Info” 
panel, “Source Info” panel, mode selection panel, and buttons panel (Figure 3.3). The 
first four panels load the main input data file, and the “Source Info” panel loads the sound 
source input file. In the “Output Data Files” panel, the output file name and log file name 
can be given. The text fields in the “Depth Info” and “Range Info” panels decide the 
depth and range information for the output data. The text fields in the “Basic Info” panel 
provide parameters for running the model, for example, FFT size. If the value of 
‘Number of Depth Points” is zero, depth grid size is determined according to accurate or 
efficient mode. If it is not zero, it is assigned to the variable, deciding the FFT size. The 
value of “Absolute Max Depth” is the maximum depth of the deep bottom. If this value 


is Zero, it is reset to the two times the maximum depth of the water/bottom interface. 


The radio buttons in the mode selection panel for depth and range are used in 
deciding the grid size of depth and range. The buttons panel has five buttons: “Run” 
button, “Plot” button, “Input File Info” button, “Clear” button, and “Help...” button. The 
“Run” button starts executing the model. After the model computations are completed, 
the output binary file is created and the “Graphic Control” frame is automatically opened. 
The “Plot” button renders to open the “Graphic Control” frame with the output binary file 
given in the text field of the output file name. The “Input File Info” button renders to 
open the “Input File Editor” frame with the input files recorded in the main input file. 
The “Clear” button clears all text fields. The “Help...” button shows the explanation 


about this frame. 


3. Input Data Editor Frame 


This frame loads five input data files, specifically sound speed profile (SSP) data 
file, bottom bathymetry data file, bottom property data file, deep bottom bathymetry data 


file, and deep bottom property data file. Like the main input file and the sound source 


ANG 


input file, these files have the same form as the ones used in the MMPE model. The 
format of the data files is described in Smith and Tappert (1993). The “Input Data 


Editor” frame has a tabbed pane with five panels as shown in Figure 3.4. Each panel has 


. 
* 


a file name text field, “Save” button, “...” button which renders a file chooser to open an 


input data file, and a text area. 
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Figure 3.3 — Model Control Frame 


User evaluation found that the “Input File Editor” frame does not work well as an 
ASCI file editor. JTextArea, a Java class in the Java standard library, provides an 
editing capability for the input files. However using this class is not desirable because it 


sets the size of rows and columns with two dimensional character arrays. Therefore, 


future work needs to avoid using JTextArea class. JEditorPane might well be used 


instead. The problems in this frame are left for future work. 
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Figure 3.4 — Input File Editor Frame 


4. Graphic Control Frame 


As described in Chapter II, in the post-processing of MMPE, there are six kinds of 
options for graphic manipulation. In WinMMPE, four options have been worked so far: 
“Single Radial,” “Single Depth,” “Single Interface” and “Source Field.” The “Graphic 
Control” frame has a panel and a tabbed pane and is shown in Figure 3.5. The output 
binary file is chosen in the panel. With the “...” chooser button the output binary file that 
users want to use can be selected and loaded. The tabbed pane has the four panels with 
each graphic option. The graphic output from each graphic option and the components of 


each panel in details are described later in this chapter. 
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Figure 3.5 — Graphic Control Frame 


5. Graphic Output Frame 


The above three frames, “Model Control” frame, “Input File Editor” frame, and 
“Graphic Control” frame are displayed inside the main frame, while the “Graphic 
Output” frame is displayed independently on the main frame and is shown in Figure 3.6. 
The utle bar of the frame shows the chosen graphic option from the “Graphic Control” 
frame and the output file name. There are four menus: “File,” “Edit,” “Window,” and 
“Help.” As shown in Figure 3.6 (a), the “File” menu has seven menu items: “New,” 
“Open,” “Save,” “Save As,” “Print,” “Page Setup,” and “Close.” Of these menu items, 
“Print,” “Page Setup,” and “Close” do work and the others do not work. In addition, 
“Edit” and “Window” menus do not work, either. They are left for future work. The 
“Help” menu shows the user’s choice from the option in the “Graphic Control” frame. It 


is shown in Figure 3.6 (b). 


Two kinds of graphic output can be launched in this frame. One is a color map 
image of a transmission loss field (Figure 3.6 (a)). If the user would like to know the 
value of a transmission loss at a point inside the image, the user can put the cursor at that 
position and then click. The value is shown at the bottom label of the graphic frame. The 
bar of color index with respect to transmission loss is also given. Another type of plot 1s 
a single line of the transmission loss (Figure 3.6 (b)). The number and the value of the 
scale in the x-axis and y-axis of the plot are decided by the value depending on the 


difference of minimum value and maximum value. 


To develop the algorithm for the color map shown in Figure 3.6 (a), a lot of 
experimental trial and error regarding treatment of RGB values was performed. The final 
color method is shown in Figure 3.7. The red, green, and blue values of the color map 
start from 0.5, 0, and 0, respectively at the minkmum TL. They finish at 0, 0, 0.5 at the 


maximum TL. 
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Figure 3.6 — Graphic Output Frame 
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Figure 3.7 — Algorithm of Color Map 


C. IMPLEMENTATION 


As described in the previous section, all of the functionality in WinMMPE is 
initiated by clicking a button in the main frame. Model execution is initiated by clicking 
the “Run” button in the “Model Control” frame. The graphic manipulation is initiated by 
the “Plot” button in the “Graphic Control” frame. For graphic image generation, each of 
four graphic options is described in this section. This section does not describe what 


kinds of data used, but rather how to manipulate the graphic outputs. 


1. Model Execution 


As soon as the model starts running, processing outputs appear in the DOS 
command window as shown in Figure 3.8. The program then produces a binary output 
file. The file has the same structure as the one from the MMPE model. Ordinarily the 
binary serialization file in Java uses big-endian floating-point representations as a default 


type in the standard Java library supported by Sun Microsystems Corp. However, in the 
Sy 





WinMMPE, the little-endian representation is used for the binary output file to be 


compatible with the Windows version of Matlab for post-processing of the MMPE model 


data. The little endian java class is referenced at http://mindprod.com/products.html - 


LEDATASTREAM. 
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2. Post-processing and Graphic Manipulation 


The “Single Radial” post-processing option in the “Graphic Control” frame is 
shown in Figure 3.5. The graphic output is produced by Eq. (2.53). It has a “Radial” 
section, a “Frequency” section, and four buttons: “Plot,” “WRML,” “Clear,” and “Help.” 
Each section has an input text field. Whenever the output file is loaded, 0.0 degrees and 
the center frequency are set as defaults in the “Radial” and “Frequency” fields, 
respectively. The “Plot” button renders a “Single Radial” graphic output. The title bar of 
the frame shows the phrase, “Single Radial,” and the output file name. The acoustic 
transmission loss field at a single radial along the depth and the range is shown in Figure 
3.9. The single line is the bottom bathymetry as a function of range. It is possible to 


show the deep bottom bathymetry as well, if it exists in the data. 
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Figure 3.9 — Single Radial Graphic Output 


The “VRML” button is only in the “Single Radial’ option. This button renders to 
create four ASCII data files for showing a VRML scene graph. They are “point” for 
point data, “color’ for color data, and “indexPoint” for coordIndex and colorIndex in 
VRML. These files are combined in a VRML file, “mmpeVRML.wrl.” The WinMMPE 
then renders a web browser set as default on the computer. The web browser must be 
able to launch the “VRML” file with the extension “.wrl.”” Then the 3D scene of the 
acoustic transmission loss field is simulated, as depicted in Figure 3.10. In fact, the data 
can be viewed with the true aspect ratio in terms of real range and real depth. Note that 
the scene of Figure 3.10 is minimized by 1/1000 for range and by 1/100 for depth. 
Currently the 3D scene does not have much meaning since it is simply a single 2D 
acoustic field in 3D space. Future work expanding the 3D model is discussed in Chapter 
V. 
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The GUI of the next option, “Single Depth”, 1s shown in Figure 3.11. The 
graphic output is produced by Eq. (2.53). It has a “Depth” section, a “Frequency” 
section, and three buttons. When the output file 1s loaded, the source depth and center 
frequency are set as defaults in the “Depth” and “Frequency” fields, respectively. The 
“Plot” button renders a graphic output of “Single Depth.” Figure 3.12 shows two cases 
of graphic outputs. Figure 3.12 (a) shows a multi-radial acoustic field and Figure 3.12 (b) 
shows a single radial acoustic field. Note that if multiple radials exist in the data file, the 


plot automatically displays the TL field of radial versus range. 


a eta eer agree 


= MMPE 2D Graphic Control Frame 20000 








j On ee ee STI en pe NT Lon ON 
| {Single Radial “Single Depth [Single Interface ‘Source Field 


————— 
Enter Depth Gn metezs positive dowravand) (50.0 | 





“Frequency 


Enter Frequency 190.0 Hz to 108.9 Hz 1 00.0 | 


———_. et - ‘ —— 
Piot Clear | Help... | 





Figure 3.11 —- The Option “Single Depth” in Graphic Control Frame 
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The GUI of the “Single Interface” is shown in Figure 3.13. The graphic output is 
produced by Eq. (2.53). It has a “Depth” section, a “Frequency” section, and three 
buttons. Unlike the “Depth” section in the “Single Depth” graphic option, it has two 
radio buttons for choosing one of either “Water/Bottom Interface” or “Bottom/Basement 
Interface,” instead of an input text field. When an output file is loaded, “Water/Bottom”’ 
and center frequency are set as defaults in the “Depth” and “Frequency” fields, 


3 


respectively. The “Plot” button renders a graphic output of “Single Interface.” Figure 
3.14 (a) shows a mult-radial acoustic field at “Water/Bottom Interface” and Figure 3.14. 
(b) shows a single radial acoustic field at “Bottom/Basement Interface,” as shown in the 


title bar of each figure. 
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Figure 3.13 - The Option “Single Interface” in Graphic Control Frame 


(a) 
(6) 





Tv 
we 
F ts La) Co) 
> © > a | ~~ 
i> fav] ww - vw 
CR ee eat eataateaea nt fale aon" mamamemamcamamamamma eeeeeteeeched 5 
' 7 omnganae nee e ‘ 
& “g : gp carn ® DOPeEHO NOE eg SOON 
? : Ne ooo paren ca evden sk me, - 
j nweees oe leamaaiel 
4 ; Fag tgs ee 
i Pyiterey ete 
<Neene Ri me nase 9 ere crzeisons mecreemaconacnte® At eneardhibesesrenenentn ticeltoetads 
ae Nit odin dad 


ep OG ed 
CMa ne a on sw Geer eres 
PQ 44 He Be eg 
rr enn ee tow 
BOON aaaagose 
PE ee nwene were «, 


ang 8 NPE spect Sah ES 


x 


eternal 
ee 


< '§ 3 ee 
Sy mes a Peteraye >» 
4 


Pe: 
ee we Seal & dytet4% bat 
eat ph ap ental Se tt 
MB ETT wan 
ae Ore rerenrees o¥ 
an ttn Oe ROI ae 
aod : 
Bee 
bh nccitintinntadin re TT 
§ men na a pinche 6 
Yr? 


: 
a> 
x: ee 
nn aR smear tty tt: 
<<“ ~ 


0 


SEE 


v 
PSial Gey 
TL Information (Click iriside thre knage} 





$ (SE 7 Ir) 


NPP ST EIT SRI ee 
ee tad x 


vp arma Sete 


= 
. 


t 
wee maLsorres SE ee wre ene 
a coatinined af Lae ‘ 
F ie ecauigagante 
wens 2eeF hr s — perme SOMO ALR My on ae 
Pertret  wemsenaane 
ane ty * oa tt oe 
maine OS" 


a 


TRANSMISSION LO. 


POR we PR ee ee 





a aie = ae a 
DES AOE din Ne SORES RENE RELA IH LAPT DOG : 


~~ ~~ 
psy a 
B 
ba aed . 
OO er | Vane RTE CP NALAY Bovrewenorevererdweres att eadanee me be aeeumeou Need Rice. ’ 


te> a t+, ™ thr 
a tt es + a) ay fw. 0% est ve 


QU4) TORINO dt 22 DE UL 





Fie Eat Wino Help 


Fie Eck Viitsiow Help 
4 
r 
7 





(b) 


+ 
> 


So 


4GE vem: 


==] 


WS 


~~, 


14 — Single Interface Graphic Output: (a2) Multi-Radial Acoustic Field 
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Single-Radial Acoustic Field 


igure 


Figure 3.15 shows the GUI of the “Source Field.” The graphic output is produced 
by Eq. (2.55). It has only “Plot” and “Help...” buttons. The “Plot” button renders the 
two kinds of graphic outputs, one of which is broadband and the other is for a single 
frequency data file. The former is shown in Figure 3.16 (a) and the latter is shown in 
Figure 3.16 (b). Note that if multiple frequencies exist in the data, this option 


automatically creates a frequency versus depth TL field plot. 
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Figure 3.15 — The Option “Source Field” in Graphic Control Frame 


40 






peg st wetieny | weld (Mult FreqPeout.bin) — a ; i ae ff «| 





TRANSMISSION LOSS (dB re 1m) 


$0 
100 
150 


(a) 


200 


DEPTH (rrp 


750 


300 


350 
100 





80 90 100 116 120 
FEQUENTY (Ha} 


TL information (Cick inside the image) 





fn Starting Fiehd (siniple freaut Bin) =. we F | | x 


File Edt Window Hey ae 
$0 
100 
150 


(b) 


200 


DEPTH (ry 


240 


300 





su 


10 70 30 49 50 40 70 86 30 100 
TL (dB re Im) 


Figure 3.16 — Source Field Graphic Output: (a) Broadband Acoustic Field; (b) Single- 
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IV. OUTPUT RESULTS 


A. INTRODUCTION 


In Chapters II and III, the implementation and data processing of the MMPE and 
WinMMPE models have been described. In this chapter, a comparison of the model 
output data, the values of transmission loss, and graphic outputs between MMPE and 
WinMMPE are made in order to confirm the accuracy of the new model. The efficiency 
of WinMMPE compared with MMPE is also examined. Specifically, this comparison is 


based upon the amount of memory used and the computational run-times. 


As examples of typical cases, two sets of environmental input data files have been 
defined. The first case is based on a simple, range-independent Pekeris waveguide. The 
water column is isospeed with sound speed 1500 m/s. The bottom is at constant depth 
250 m, and consists of a homogeneous, fluid half space with sound speed 1600 m/s, 
density 1.2 g/cm’, and attenuation 0.15 dB/km/Hz. The source is at a depth of 100 m and 
transmits a CW tone of 100 Hz. The second case is based on real environmental 
measurements taken during the ONR-sponsored Primer experiment off the east coast of 
New Jersey in 1996 (Lynch et al, 1997). The source was at a depth of 270 m and the 
acoustic energy propagated up the continental slope onto the shelf. The acoustic 
frequency for this run was 400 Hz. Hereafter, the first test case (Pekeris waveguide) will 


be referred to as Case I, and the second test case (Primer) will be referred to as Case I. 


As shown in Table 4.1, Case I has only one profile within each of the 
environmental input data files. The total number of environmental data values are 20. By 
contrast, Case II has 319,872 different sound speed profiles and 29,952 different 


bathymetry data points. The total number of environmental data values are 699664. 
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a = 


128 radials X 

Sound Speed Profile Data (2 data values) I 49 ranges X 
51 depths 

128 radials X 


Bottom Property Data (7 data values) St 
ee 
Deep Bottom Property Data (7 data values) 


Table 4.1 — The Number of Input Data Profiles 







Both the input data and specification of parameters within “pefiles.inp” decide the 
matrix of the output data from the model. Specifically, all radials defined by the 
environment and all frequencies defined in the source input file are output. However, 
only the number of depth and range points requested between certain limits defined 
within “pefiles.inp” are stored. Thus, Case I has a two-dimensional matrix, as shown in 


Figure 4.1 (a), of size 103x200. Case I, on the other hand, has a three-dimensional data 


matrix, as shown in Figure 4.1 (b), of size 128*128x200. 


103 depthsteps (a) 


200 Range Steps 





128 Radial Steps 


128 Depth Steps (b) 


Figure 4.1 — Matrix dimensions of output data for (a) Case I and (b) Case II. 
44 


B. COMPARISON OF OUTPUT DATA 


The Fortran version of MMPE primarily uses single precision floating point 
values with each data occupying 4 bytes. In contrast, the Java implementation of 
WinMMPE uses 8-byte, double precision data types. There are two reasons for using 
double precision. One is that the default variable type in Java is double precision. If a 
developer would like to use single precision variables instead, the cast operator with the 
“float” command must be added in the definition of the value. Another reason to use 
double precision is that more precise outputs can be obtained. However, the output data 


from both models is single precision. 


This section shows the comparison of output data from the MMPE model and the 
WinMMPE for the given environmental input data. Some of the values to be examined 
will be the maximum difference of the PE field function and the maximum difference of 
transmission loss (TL) between MMPE and WinMMPE for each graphic. In addition, 


more general comparisons of the graphic outputs from each model are made. 


1. Single Radial 


We begin by examining the differences between single radial results. The PE 
field function as defined in Eq. (2.16) is a complex valued vector. As shown in Table 
4.2, the maximum amplitude difference of the PE field function for Case I is 
2.4780x10° while the maximum difference in TL is 0.0054 dB re lm. These 
differences are negligible and can be attributed to the differences between the single 
precision and double precision implementations. Figure 4.2 shows the “Single Radial” 


graphic outputs from MMPE using the Matlab post-processing routines and WinMMPE. 


In Case II, the amplitude of the maximum difference of the PE field function is 


1.4950x10~ while the maximum difference in TL is 7.5633 dB re 1m (Table 4.2). This 
large TL difference suggests there may be problems between the solutions. However, 
upon closer inspection, it is found that all of the significant differences occur in regions 


where the PE field function is extremely small, i.e. at discrete locations within the bottom 
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volume. Table 4.3 lists all such points in range and depth where the TL difference 
exceeds 3 dB. In total, there exist 8 points over 128 depth and 200 range points. Figure 
4.3 shows the “Single Radial” graphic outputs from MMPE and WinMMPE along the 
central radial. These plots confirm that the locations of the discrepancies are in regions 
of extremely low field strength. Again, this could be due simply to the difference 
between single and double precision implementations. Overall, these plots indicate that 


WinMMPE is working properly and producing accurate results. 


Maximum Difference 
( [MMPE - WinMMPE| ) 
PE ~ 
— Field Function 2.4780x 10 
Case I = 
(dB re 1m) 0.0054 
PE 


—4 
Field Function 1.4950 10 


TL 
Chia 7.5033 


Table 4.2 — Maximum difference between the MMPE Model and the WinMMPE Model 
for Single Radial. 






Case II 









Table 4.3 — Positions where TL difference exceeds 3 dB for Case II. 










46 


r "3 f — “ 
= ah De 
= - ——=s i _ - 


—s i 


— 
ge Hy eens meee kr este ey oy on ae 


(a) 





TRAMCAISSOY LOB (8 re tim 


(b) 


CEPTH iy 





SAME sbrrit 


TL Indormation (Gick inside the image) 





Figure 4.2 — Single Radial graphic output for Case I: (a) from MMPE; (b) from 
WinMMPE. 
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Figure 4.3 - Single Radial graphic outputs for Case HH: (a) from MMPE; (b) from 
WinMMPE. 
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2. Single Depth 

We now consider the data produced by selecting the single depth output option. 
The maximum amplitude difference of the PE field function for Case I is 2.1393x10~* 
while the maximum difference of TL is 0.0010 dB re 1m. These values are recorded in 
Table 4.4. Again, these differences may be considered negligible. Figure 4.4 shows the 
“Single Depth” graphic outputs from MMPE and WinMMPE at the source depth. 


For Case II, the maximum amplitude difference of the PE field function is 0.0016 
while for TL the difference is 6.5112 dB re 1m. As before, the positions in space where 
the TL difference exceeds 3 dB are determined and listed in Table 4.5. In this case, there 
exist 9 points total over 128 radial and 200 range points. Figure 4.5 shows the “Single 
Depth” graphic outputs from MMPE and WinMMPE at the source depth. As expected, 


these large differences only occur where the acoustic field has extremely low values. The 


similarity of the plots indicates that WinMMPE is providing accurate solutions. 
Test Case Dan Maximum Difference 
(|MMPE - WinMMPE] ) 
re zr 
Field Function 2.139310 
Case | 


ARE. 
(se We 0.0010 
RE 


Field Function 0.0016 


ay 
eas 6.5112 


Table 4.4 — Maximum difference between MMPE and WinMMPE for Single Depth. 


ee ee 2 


Case II 





Table 4.5 — Positions where TL difference exceeds 3 dB for Case II. 
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Figure 4.5 - Single Depth graphic outputs for Case II: (a) from MMPE; (b) from 
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3. Single Interface 


As was described in Chapter IJ, there are two radio buttons in the “Single 
Interface” graphic frame which allow the user to extract data either from the 
“Water/Bottom” interface or the “Bottom/Basement’” interface. Because the 
environments examined here do not contain the deep bottom/basement interface, we will 
only examine the “Water/Bottom” interface option in this section. Table 4.6 shows that 
the maximum amplitude difference of the PE field function for Case I is 1.545110 
while the difference in TL is 3.6884x10~ dB re Im, both of which may be considered 
negligible. Figure 4.6 shows the “Single Interface” graphic outputs from MMPE and 
WinMMPE at the “Water/Bottom” interface. 


For Case II, the maximum amplitude difference of the PE field function is 0.0017 
and the maximum difference in TL is 0.5484 dB re 1m, as stated in Table 4.6. Because 
we are explicitly examining values along the interface, we do not encounter the extremely 
small values of the acoustic field that were problematic before. Furthermore, the largest 
differences occur in the TL “nulls”, which may also be affected by the double precision 
implementation. As a result, even these differences for Case II may be considered 
negligible. Figure 4.7 shows the “Single Interface” graphic outputs from MMPE and 
WinMMPE at the “Water/Bottom” interface. 


Maximum Difference 
(|MMPE ~ WinMMPE| ) 
Acoustic = 
Field Function 1.5451X10 
Case I = 
(dB re 1m) 3.6884xX10~ 


Acoustic 
Field Function 0.0017 


iIE 


Table 4.6 ~ Maximum difference between MMPE and WinMMPE for Water/Bottom 
Interface. 
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Figure 4.7 — Water/Bottom Interface graphic outputs for Case II: (a) from MMPE; (b) 
from WinMMPE. 
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4. Source Field 


We now examine the differences in the computed starting fields at zero range. 
Table 4.7 shows that the maximum amplitude difference of the PE field function for Case 
Lis 2.5333x107 and the maximum difference in TL is 0.0664 dB re 1m. For Case II, the 
maximum amplitude difference of the PE field function is 3.1593x10~° while the 
maximum difference in TL is 0.2548 dB re Im. In both cases, the differences may be 


considered negligible. 


Figure 4.8 displays the “Source Field” graphic outputs from MMPE and 
WinMMPE at the initial range for Case I. The corresponding graphic outputs for Case II 


are provided in Figure 4.9. 


Test Case Dara Maximum Difference 
( IMMPE — WinMMPE] ) 
| Acoustic = 
Field Function 2.5333X10 
Case I 


THE 
(dB re lm) ge 


Acoustic 
Field Function 
TL 
EB 0.2548 
Table 4.7 — Maximum difference between MMPE and WinMMPE for Source Field. 


3.1593x10° 
Case II 
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Figure 4.8 — Source Field graphic outputs for Case I: (a) from MMPE; (b) from 


WinMMPE. 
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Figure 4.9 — Source Field graphic outputs for Case II: (a) from MMPE; (b) from 
WinMMPE. 


= 


C. EFFICIENCY 


The environment defining Case I was very simple with only a small number of 
input parameters. It was also range-independent, so the propagation operator only needed 
to be defined once and could be used throughout the calculation. Case II, on the other 
hand, had a very large number of environmental input parameters and was range- 
dependent. This required an update of the propagation operator at every range step. 
Therefore, it is to be expected that Case II will take longer to mun than Case I. 


Furthermore, the amount of memory used for Case II should be larger than for Case I. 


This section describes the efficiency of the MMPE Model and the WinMMPE 
Model as both a propagation model and a graphic tool in terms of main memory usage 
and running time. All testing was run on a single system with the following 


specifications: 

e Intel Pentium II Processor, 366 MHz; 

e Main Memory, 392 Mbytes; 

e Operating System, Microsoft Windows 2000; 

e Digital Visual Fortran Compiler, Standard Edition, Version 5.0.A; 
e Matlab, Version 5.3.1 (R11.1) with Signal Processing Toolbox; 

e Java2 SDK, Standard Edition, Version 1.3.0; 

e Forte For Java, Community Edition, Version 1.0, Update Release 1. 


The main memory on this system is adequate to avoid the need for swapping with the 
hard disk. Thus, although faster systems are available today, this should provide a 


reasonable platform to compare model performance. 
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1. Main Memory 


As mentioned in Chapter I, programmers using standard Fortran must define the 
size of data arrays at the beginning of the program. This fixes the size of all arrays 
throughout the calculation, regardless of their usage. More recent versions of Fortran 
compilers can handle dynamic memory allocation, but this will not be available on all 
platforms. In this analysis, therefore, we will assume the Fortran program has a fixed 
size of data arrays for input data. If users want to use larger input data sets, the program 
must be re-coded to accommodate the larger data array sizes. In general, when a very 
simple data set is used which only occupies a fraction of the fixed allocated space, the 


rest of the memory is wasted. 


In Java, however, programmers create the data arrays dynamically. This allows 
the program to allocate memory space only when it needs it and only as much as it needs. 
This provides for a much more efficient use of system resources than the Fortran version. 
Even though the Java data are 8-byte double precision variables as compared to the 
mostly 4-byte single precision variables in the Fortran version, this dynamic allocation of 


memory is expected to significantly reduce the memory requirements of the program. 


The amount of memory used in each version of the model is shown in Table 4.8. 
These values were obtained from the “Task Manager” program of the operating system. 
Although these numbers do not represent exact values, they do provide a reasonable 
estimate of the memory usage of each model. For MMPE, both Case I and Case I 
roughly took 120 Mbytes. The memory usage of WinMMPE is due to both basic model 
functions (e.g., the GUI objects), and the dynamic allocation for run-time variables. The 
basic functions took about 18 Mbytes. The dynamic allocation portion took about 4 


Mbytes for Case I and 17 Mbytes for Case II. 


3E, 












Model 
MMPE (Fortran) = 120 Mbytes = 120 Mbytes 


= 18 Mbytes (basic) = 18 Mbytes (basic) 
+ 4 Mbytes (dynamic) + 17 Mbytes (dynamic) 





WinMMPE (Java) 






i 





Table 4.8 — Main memory used in the models for Case J and Case II. 


For the graphical outputs and post-processing, the amount of memory used in 
WinMMPE is almost the same as that used in Matlab. This is because Matlab uses 
double precision and dynamically declares the size of data arrays with given sizes from 


the output binary file. 


2. Running Time 


Because the model is based on the split-step Fourier algorithm, a significant 
fraction of the run-time is due to the fast Fourier transforms (FFT) between depth and 
vertical wavenumber. In Case I, the number of depth points, 1.e., number of terms in the 
FFT, is 1024. In Case II, the number of depth points is 2048. For each range step, two 
FFT’s must be performed. Case I requires 222 range steps as compared to 15490 range 
steps for Case II]. Furthermore, for multiple radials or multiple frequencies, the same 
number of FFI’s must be performed for each range calculation. The total number of 
FFT’s for some calculations can then become significant. Table 4.9 shows how many 
times the FFT module in the model is called for each case. Note that Case II has over 


8900 times more FFT calls than Case I. 


a 

ee re a 

ranges 
fe Tin 


| Total = Radials x Frequencies x Ranges x 2 3965440 
Table 4.9 ~ Number of FFT calls in the PE/SSF algorithm for Case I and Case IL. 
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Table 4.10 shows the run-times of the models for Case I and Case II. This 
measurement was achieved by watching the timer in the operating system. In some 
cases, it is hard to tell how long the model takes to complete, especially in the case of 
small input data and relatively few range steps. For Case I, the MMPE Model spent less 
than a second and the WinMMPE Model spent less than two seconds. Thus, it is hard to 
compare the true relative run-times. Case II, however, took a much longer time to run. 
The MMPE Model took about 3 hours and 30 minutes. In contrast, the WinMMPE 
Model took 6 hours and 30 minutes. It appears, therefore, that WinMMPE takes about 
twice as long to complete a calculation as MMPE. This factor of two can be a significant 
issue for very large calculations. It is quite possible that the problem is, in part, due to 
the Java code method of calling the “getFFT( )” in the “FFT.java” class in order to 
implement the FFT. The Java class was designed like the Fortran FFT subroutine in the 
MMPE Model, which is based on a standard numerical technique (Press, et al., 1992). It 
is recommended that the FFT routine be re-coded in an object-oriented programming 


style. This will be discussed further in Chapter V. 


Model 


— 


Table 4.10 — Run-times of the MMPE and WinMMPE Models for Cases I and I. 





The run-time of the post-processing and graphic presentation is listed in Table 
4.11. For Case I, the Matlab routines spent a maximum of 2 seconds generating the 
‘Single Radial” output while the WinMMPE Model spent only 3 seconds. The run-times 
for Case II are nearly the same for the two models. For the “Single Depth” graphic 
output, the run-time of the WinMMPE Model is typically 4 seconds. The Matlab 
routines, on the other hand, consistently took about 41 seconds to complete. For the 
“Single Interface” graphic output, WinMMPE generally took about 24 seconds while 
Matlab still took about 41 seconds to complete. The run-time of the other optional 


graphic output, “Source Field,” was not significant. 
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It should be noted that these run-times do not account for the time taken by the 
user to run the propagation model and the post-processing routines separately. Time is 
also taken to manually enter parameter values on the command line interfaces for the 
MMPE Model and the associated Matlab routines, as well as the time taken to fill in the 
parameter fields within the WinMMPE Model. The WinMMPE Model at least has the 
advantage of a unified program within which GUI’s exist to control the propagation 
model, post-processing, and graphic presentations. As a result, the operator’s skill and 


familiarity with the programs can affect the run-times. 














Graphic Tool 


WinMMPE (Java) 


Table 4.11 — Run-times of the post-processing and graphic presentations for Case I and 
Case II 













Single Radial: = 3sec 
Single Depth: = 41sec 
Single Interface: = 45sec 
Source Field: < lsec 


Single Radial: < 2sec 
Single Depth: < lsec 
Single Interface: < lsec 
Source Field: < lsec 







Single Radial: = 3sec 
Single Depth: = 4secs 
Single Interface: = 24secs 
Source Field: < lsec 


Single Radial: < 3sec 
Single Depth: < Isec 
Single Interface: < lsec 
Sourcedsieldau< lisec 
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V. CONCLUSIONS & RECOMMENDATIONS 


There are a variety of underwater acoustic propagation models available for the 
general user these days. These include models based on ray theory, normal modes, 
wavenumber integration techniques, and parabolic equation models. The ONR- 
sponsored web site, http://oalib.saic.com, provides copies of many current research 
models used in the professional community. However, these models are not usually very 
user-friendly for the untrained. None of the models examined provide a convenient 
eraphical user interface and few have been developed in an object-onented language. To 
address these issues and provide a more user-fnendly operating environment, a Windows 


version of the MMPE Model, referred to as the WinMMPE Model, was developed. 


The MMPE Model is really composed of two parts. The propagation component 
was written in Fortran, while the post-processing and graphical presentations were treated 
using Matlab routines. The WinMMPE Model, on the other hand, incorporated both 
parts into a single program written in Java. In order to compare the accuracy and 
efficiency of the WinMMPE Model with the previous MMPE Model, two test cases were 
defined. The first case was a simple, range-independent, single radial Pekeris waveguide 
with very few environmental input parameters. The second case was a much more 
complicated, range-dependent, 3D environment with multiple radials based on measured 


environmental data taken near the Mid-Atlantic Bight off the coast of New Jersey. 


A comparison of the data indicated that WinMMPE was producing results nearly 
identical to MMPE. The largest differences were observed in regions where the acoustic 
field strength was extremely small and so could be due to simple round-off error 
differences. This is certainly likely since the default data type in the Java implementation 
of WinMMPE is double precision while the default type in the Fortran implementation of 
MMPE is only single precision. In regions where the acoustic field strength was large, 


the differences were only fractions of a dB. 
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A comparison of the efficiency of each model in terms of the ieiiow usage and 
run-time was also presented. Due to the initialization of all variable array sizes, the 
Fortran-based MMPE Model used the same amount of memory for both test cases. In 
contrast, the Java-based WinMMPE Model dynamically allocated the necessary memory. 
The MMPE Model always required about 120 Mbytes, as compared to the roughly 25 - 
30 Mbytes needed to run the propagation component of the WinMMPE Model. With 
respect to the run-times, is seemed that the WinMMPE Model generally took about twice 
as long to compute the propagation solution as the MMPE Model. This could be due in 
part to the increase from single precision to double precision. It was also suggested that 
the FFT routine was not well suited for implementation within Java. Finally, it was 
found that there was essentially no difference in run-time for the post-processing and 


graphic presentation of results. 


There are five essential recommendations for future work. The first 
recommendation concerns the FFT routine. As suggested above, the FFI Java class 
should be re-created or replaced with an existing FFT Java open source on the Internet. 
There is currently no FFI Java class in the standard library provided by Sun 
Microsystems Corp. However, more and more class libraries are becoming available 
over the Internet. One of them is_ provided on the web site 
http://www.fusion.kth.se/courses/jbone. Complex classes must be combined with it as 
we.. Therefore, some parts of the algorithm must be changed since the treatment of real 


and imaginary data must be combined to treat a complex variable. 


Second, there are additional graphic output modules that need to be added in 
WinMMPE. The Matlab routines of the MMPE Model have two more options than the 
WinMMPE Model, specifically “Single Range” and “Travel Time”. The first of these 
provides another view of the data to help the user evaluate the acoustical environment. 
The latter option is of great interest for any future development since it provides for the 
Founer synthesis of broadband data to produce predictions of pulse propagation in the 
time domain. There are also options that allow the user to produce plane-wave 


beamformed results, thereby providing the user with arrival angle versus time data. 
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Furthermore, there are some menu items in the graphic frame that do not yet work, as 


mentioned in Chapter III. 


Third, the problems associated with the “Input File Editor” frame were described 
in the “Design and Functions” section of Chapter IJ. Even though that frame was not a 
major work in this thesis, it is recommended that the frame be changed to work like the 
well-functioned standard Windows editor. Additionally, the functions of the menu items 


in the “Edit” menu of the main frame should be upgraded to work as expected. 


Fourth, it would be a great service to the general acoustics community if the 
WinMMPE Model was available over the Internet using Java applets. The WinMMPE 
Model might provide useful information about the underwater acoustic environment for 
scientists who are not necessarily acoustics experts, or even for the general public 
interested in scientific modeling and visualization. It could even be used by schools to 


interest school kids in science and the ocean environment. 


Finally, 3D visualization of underwater sound propagation can be an interesting 
and useful project in the future. Holliday (1998) worked on real-time 3D sonar modeling 
and visualization with a ray-based model. His work provided an application- 
programming interface (API) for real-time 3D sonar modeling and visualization of sound 
propagation. Future work could consider a similar 3D sonar visualization with the 
WinMMPE Model. Figure 3.10 was used to illustrate the scene of an acoustic field for a 
single radial. This simple example was not, however, a true 3D acoustic field since it 
only provided a single plane of 3D structures to visualize. Hopefully, Figure 3.10 can be 
a Starting point to develop the 3D visualization of the actual 3D acoustic data provided by 


WinMMPE. 
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